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Abstract: Assuming 1^(10890) as a P-wave 6g-scalar-diquark 6g-scalar-antidiquark o 
tetraquark state, the mass of Yj>( 10890) is computed in the framework of QCD sum rule 
method. Technically, contributions of operators up to dimension six are included in the 
operator product expansion (OPE). The numerical result 10.88 ± 0.13 GeV for 1^(10890) 
agrees well with the experimental value, which favors the P-wave [ftgjffeg] tetraquark con- 
figuration for 1^(10890). In the same picture, the mass of y(4360) is calculated and the 
result 4.32 ±0.20 GeV is compatible with the experimental data, which supports Y(4360)'s 
P-wave [eg] [eg] structure. 
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1. Introduction 

The observations of T(1S)tt + tt^ and Y(2S)ir + ir~ states near the T(55) resonance [|l], |2[ 

have attracted great theoretical attention |J. However, there are still some puzzles on 

the anomalously large rates and the way to describe distribution shapes and the helicity 

angle. Recently, Ali et al. {§ identify V fe (10890) with the state Y [bq] (10900) [§ and interpret 

Yb(10890) as a P-wave tetraquark state. In this way, a dynamical model for decays 

Yf, — > Y(lS')7r + 7r~, T(2S)tt + tt~ is presented, which provides excellent fits for the decay 

distributions. Therefore, it is interesting to investigate whether 1^,(10890) could be a 

tetraquark state. Undoubtedly, the quantitative description of Yt(10890)'s properties like 

mass is helpful for understanding its structure, but it is difficult to extract the hadronic o 

spectrum information from the simple QCD Lagrangian. That is because low energy QCD 

involves a regime where it is futile to attempt perturbative calculations and one has to treat 

a genuinely strong field in nonperturbative methods. However, one can apply QCD sum 

rules H (for reviews see [[/], ^, |9], and references therein), which are a nonperturbative 

formulation firmly rooted in QCD. From the above reasons, we devote to study 1^(10890) 

with QCD sum rules in this work. 

Additionally, BABAR Collaboration observed a broad structure Y(4325) in the process 
e+e" -> 7/ Si? 7r + 7r-'0(2s) at 4324 ± 24 MeV with a width 172 ± 33 MeV Q. Latterly, 
Belle Collaboration reported the charmoniumlike state Y(4360) in e + e~ —> ir + ir~ip(2s) at 



4361 ±9±9 MeV with a width of 74±15±10 MeV |2|. The mass of V(4325) is close to that 
of y(4360), and the main difference between them is their widths. It seems very difficult to 
observe these two structures simultaneously because of the large width of Y(4325). They 
could be the same structure and the width difference may be due to the experimental error. 



In Ref. [13], Liu et al. have tried to perform a combined fit to e + e~ — > tt + tt^^(2s) cross 
sections measured by the BABAR and Belle experiments. In this work, we assume y(4325) 
and y(4360) are exactly the same resonance for simplicity. On Y(4360), there have already 



been some theoretical works 14, 15, fig] . From QCD sum rules, Albuquerque et al. arrive 



at M = 4.49 ± 0.11 GeV, adopting the current [cq]s=o[cq]s=i + [cq]s=i[cq]s=o [0] ( for the 
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concise review on multiquark QCD sum rules, one can see |jl8[| ). At present, we would like 
to study whether Y(4360) could be a P-wave [eg] [eg] tetraquark state. 

The paper is organized as follows. In Sec. ^, the QCD sum rule for the tetraquark state 
is introduced, and both the phenomenological representation and QCD side are derived, 
followed by the numerical analysis to extract the hadronic masses in Sec. |3[ Section |] is a 
brief summary. 



2. The tetraquark state QCD sum rule 

In the tetraquark interpretation, Y\Q q ^ is a J PC = 1 diquark-antidiquark state, having 
the flavor content [QgHQg]- Its spin and orbital momentum numbers are: S'fQq] = 0, 
^IQg] = 0> S[Q q ][Qq] = 0, and -^[Q ? ][Qg] = 1 [ 19 [ - For the interpolating current, a derivative 
could be included to generate ^[Q g ][Qg] = 1- Presently, one constructs the tetraquark state 
current from diquark-antidiquark configuration of fields, while constructs the molecular 
state current from meson-meson type of fields. Although these two types of currents can 
be related to each other by Fiertz rearrangements, the relations are suppressed by color 
and Dirac factors [18]. It will have a maximum overlap for the tetraqurk state using the 
diquark-antidiquark current and the sum rule can reproduce the physical mass well. Thus, 
the following form of current could be constructed for the P-wave [Qg][Qg], 

f = e ahc e de MlClMB^q d ^CQ T e ). (2.1) 

Here the index T means matrix transposition, C is the charge conjugation matrix, 
denotes the covariant derivative, and a, 6, c, d, and e are color indices. 
The mass sum rule starts from the two-point correlator 

IT"V) = i J d?xe*' x {0\T\j»(x)j^(p)]\0). (2.2) 

Lorentz covariance implies that the correlator Q2.2| ) can be generally parameterized as 

n^(g 2 )=^-rjnW( 9 2 ) + ^n(°)( ? 2 ). (2.3) 

The part proportional to g^ v is chosen to extract the sum rule here. In phenomenology, 
nW(g 2 ) can be expressed as a dispersion integral 



IlW(g 2 ) = [ 9 J „ + - / ds— =-! — + subtractions, (2.4) 



[A^] 2 i r ^ imn( 1 )P he n (g) 

Mfj-q 2 vr J S0 s-q 2 

where Mh denotes the mass of the hadronic resonance. In the OPE side, nW(g 2 ) can be 
written in terms of a dispersion relation as 

roc n OPE M 
n^VH/ ds P - W (2.5) 

J Ami S ~ q 

where the spectral density is given by 

p 0FE (s) = -ImU^is). (2.6) 

7T 



After equating the two sides, assuming quark-hadron duality, and making a Borel trans- 
form, the sum rule can be written as 



[ A (l)]2 e -M|/M» 



dsp 0PE (s)e-°/ M2 . 



(2.7) 



To eliminate the hadronic coupling constant A^ 1 ) , one reckons the ratio of derivative of the 
sum rule and itself, and then yields 



M 



H 



4m?, 



dsp 



OPE sg - s /Af 2 



dsp OPE e -s/M* 



(2- 



To calculate the OPE side, we work at leading order in a s and consider condensates 
up to dimension six with the same techniques in Refs. |2(], 21]. To keep the heavy-quark 
mass finite, one uses the momentum-space expression for the heavy-quark propagator, 
including two and three gluons attached expressions given in Ref. p2| . The light-quark 
part of the correlator is calculated in the coordinate space and then Fourier-transformed 
to the momentum space in D dimension. The resulting light-quark part is combined 
with the heavy-quark part and dimensionally regularized. It is defined as r(rriQ,s) = 
(q + /3)mg — a[3s and K(a,/3) = 1 + a — 2a 2 + /3 + 2a/3 — 2/3 2 . The spectral density is 
written as 



p '( s ) = p P ert ( s ) + p<w>( s ) + p^)\s) 
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The integration limits are given by a m i n = (1— Jl — 4mg/s)/2, a ma;c = (1+y 1 — 4mg/s)/2, 

and /3 m j n = amg/(sa — mg). Note that the next-to-leading order corrections are not in- 
cluded here, for which one needs to consider the renormalization of the current 23]. This 
procedure is undoubtedly complicated and tedious, since the renormalization of the cur- 
rent will raise the operator-mixing problems. Especially for the multiquark system, many 
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operators will mix under renormalization. Actually, a lot of hard calculations already need 
to be done even if one works at leading order for the difficulties of tackling with the massive 
propagator diagrams. Under such a circumstance, it is expected that one could obtain a 
trusty sum rule working at leading order in a s , and it has been tested to be feasible for 
many multiquark states [18]. To improve on the accuracy of the QCD sum rule analysis, 



it is certainly meaningful to consider the next-to-leading order corrections, which may be 
included in some further work after fulfilling a burdensome task. 

3. Numerical analysis 

In this Section, the sum rule ( p. 8] ) will be numerically analyzed. The input values are taken 
as m c = (1.23 ± 0.05) GeV, m b = (4.20 ± 0.07) GeV, pi (qq) = -(0.23 ± 0.03) 3 GeV 3 , 



(gqa ■ Gq) = m 2 (qq), m 2 = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 , and (g 3 G 3 ) = 0.045 GeV 6 
H, [17], [0|. Complying with the standard criterion of sum rule analysis, the threshold so 
and Borel parameter M 2 are varied to find the optimal stability window. In the QCD sum 
rule approach, there is approximation in the OPE of the correlation function, and there is a 
very complicated and largely unknown structure of the hadronic dispersion integral in the 
phenomenological side. Therefore, the match of the two sides is not independent of M 2 . 
One expects that there exists a range of M 2 , in which the two sides have a good overlap and 
the sum rule can work well. In practice, one can analyse the convergence in the OPE side 
and the pole contribution dominance in the phenomenological side to determine the allowed 
Borel window: on one hand, the lower constraint for M 2 is obtained by the consideration 
that the perturbative contribution should be larger than the condensate contributions, to 
keep the convergence of the OPE under control and insure that one does not introduce a 
large error neglecting higher dimension terms; on the other hand, the upper limit for M 2 
is obtained by the restriction that the pole contribution should be larger than the QCD 
continuum contribution, to guarantee that the contributions from high resonance states 
and continuum states remains a small part in the phenomenological side. Meanwhile, the 
threshold parameter ^/~sq is not completely arbitrary but characterizes the beginning of the 
continuum state. On all accounts, it is expected that the two sides have a good overlap in 
the determined work window and information on the resonance can be safely extracted. 

Concretely, the comparison between the pole and continuum contributions from sum 
rule ( |2~7| ) for 1^(10890) for t/sq = 11.5 GeV is shown in FIG. 1, and its OPE convergence 
by comparing different contributions is shown in FIG. 2. In detail, the perturbative con- 
tribution versus the total OPE contribution at M 2 = 9.5 GeV 2 is nearly 63%, and the 
ratio increases with M 2 to insure that the perturbative contribution can dominate in the 
total OPE contribution when M 2 > 9.5 GeV 2 . On the other side, at M 2 = 10.5 GeV 2 the 
relative pole contribution is approximately 51%, which descends along with M 2 to guar- 
antee the pole contribution can dominate in the total contribution while M 2 < 10.5 GeV 2 . 
Thus, the regions of s and M 2 for Y b ( 10890) are taken as ^/sq = 11.4 ~ 11.6 GeV and 
M 2 = 9.5 ~ 10.5 GeV 2 . For Y(4360), the comparison between the pole and continuum 
contributions from sum rule ( |2.7D is shown in FIG. 3, and its OPE convergence by com- 
paring different contributions is shown in FIG. 4. From the similar analyzing processes, 



the regions of sq and M 2 are taken as y/so = 4.8 ~ 5.0 GeV and M 2 = 2.6 ~ 3.6 GeV 2 for 
y(4360). The corresponding Borel curves to determine masses of Y"b(10890) and Y(4360) 
from sum rule (|2.8j ) are shown in FIG. 5 and in FIG. 6, respectively. Finally, we obtain 
10.88 ±0.13 GeV for Y b (10890) and 4.32 ±0.20 GeV for Y(4360). For Y(4360), our central 
value is closer to the experimental data comparing with the prediction 4.49 ± 0.11 GeV in 
Ref. 1 17], however, the uncertainty of our result is larger. 
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Figure 1: The solid line shows the relative 
pole contribution (the pole contribution di- 
vided by the total, pole plus continuum con- 
tribution) and the dashed line shows the rel- 
ative continuum contribution for Y"&(10890). 



Figure 2: The OPE convergence is 
shown by comparing the perturbative, 
quark condensate, four-quark condensate, 
mixed condensate, two-gluon condensate, 
and three-gluon condensate contributions for 
Y b (10890). 
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Figure 3: The solid line shows the relative 
pole contribution (the pole contribution di- 
vided by the total, pole plus continuum con- 
tribution) and the dashed line shows the rel- 
ative continuum contribution for Y(4360). 




Figure 4: The OPE convergence is shown 
by comparing the perturbative, quark con- 
densate, four-quark condensate, mixed con- 
densate, two-gluon condensate, and three- 
gluon condensate contributions for y(4360). 
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Figure 5: The dependence on M 2 for the Figure 6: The dependence on M 2 for the 
mass of Y&(10890) is shown. mass of Y(4360) is shown. 

4. Summary 

In the tentative P-wave [bq] \bq] configuration, the QCD sum rule method has been employed 
to compute the mass of Yj,(10890), including the contributions of operators up to dimension 
six in the OPE. The numerical result 10.88±0.13 GeV for Y" 6 (10890) is well compatible with 
the experimental data, which favors the P-wave tetraquark configuration for Yj(10890). In 
the same picture, the mass of Y(4360) has been calculated to be 4.32 ± 0.20 GeV, and 
the result is in agreement with the experimental value, which supports its P-wave [eg] [eg] 
configuration. We expect the results could be helpful to understand the structures of these 
states. For further work, one needs to take into account other dynamical analysis to identify 
the structures of hadrons. 
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